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Abstract 

Motivated by existing and recent data on possible neutrino oscillations, we pro- 
pose a model of four light neutrinos (i/g, f^, i^t; and a singlet z^g) with a pattern of 
masses and mixing derivable from a discrete symmetry and the seesaw mechanism. 
Atmospheric neutrino oscillations occur between i/^ and Vr as pseudo-Dirac partners; 
whereas solar neutrino oscillations occur between z^e and vs, a linear combination of 
which is massless. Additional oscillations may occur between i^e and to account for 
the recent observation of the LSND (Liquid Scintillator Neutrino Detector) experiment. 



The physics of hght neutrinos appears poised on the edge of major discoveries. Many 
hints have accumulated over the past years towards nonzero masses and flavor mixing of these 
special elementary fermions without charge. They include the solar neutrino deficit []I|] , the 
atmospheric neutrino anomaly [0], the need for a cosmological hot dark matter component [0], 
and finally the excess of z/^ —>■ events observed recently by the LSND (Liquid Scintillator 
Neutrino Detector) experiment @] . If all these strands are put together in a total picture, 
along with the nonobservations of neutrinoless double /?-decay|^ and any depletion of reactor 
antineutrinos (such as from the Bugey experiment 0), one is led naturally to a variation of 
a simple scenario already proposed^]. We offer in this paper a theoretical understanding of 
this specific scenario in terms of a hierarchical seesaw model of the masses and flavor mixing 
of four light neutrinos with an assumed discrete Z5 symmetry. 

We recount the constraints from the above observations. In a two-oscillator picture in- 
volving z/g, the solar neutrino deficit favors either of two matter-enhanced oscillation solutions 
with bml^ ~ 10-*^ to 10-5 eV^, ^m^2e^^ ~ 5 x IQ-^ and with bml^ ~ 10"*^ to 10"^ eV^, 
sin^20ea > 0.4 or a vacuum oscillation solution with brv?^^ ~ IQ-^^ eV^, sin^ 2^ea > 0.75. 
Here 5mg„ is the difference of the squares of the two neutrino masses and Q^^ is the mix- 
ing angle between z/g and another light neutrino v^. The atmospheric neutrino anomaly 
can be understood in a two-oscillator scenario involving z/^ and another neutrino with 
(5m^^ ~ 10^^ eV^ and sin^26'^^ = 0(1). The LSND results, on the other hand, suggest 
(5mg^ ~ 0.5 — 10 eV^ and sin^ 2^6^* ~ (6 ± 3) x 10^^ once again in a two-oscillator picture. 
It is clear that with only the three known neutrino flavors (z/g, z^^, Vt\ we cannot explain all 
three 5m^'s. The orders of magnitude of the latter are too disparate to be explained even 
by invoking three-flavor oscillations. We are thus prompted by the need to explain all of the 
above to add a singlet neutrino z/5 which is by itself noninteracting, but will be allowed to 
mix with the other neutrinos. In this way, we are also able to consider the cosmologically 
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desirable requirement that the masses of the usual three neutrinos sum up to about 5 eV. 

With four neutrinos as suggested above, there is a simple hierarchical situation involving 
only two-oscillator scenarios!^. Let z/^ and be nearly degenerate with masses of about 
2.5 eV each and close to maximal mixing. The 5m? here is about 10~^ eV^. The singlet 
neutrino z/5 has a mass of order 10^'^ eV whereas z/g is very much lighter but mixes with z/5 
as well as z/^ by small amounts. In this scenario, because v^. is lighter than z/^, there is a 
potential conflict with rapid neutron capture (r-process) in Type II supernovaeQ. However, 
if the hierarchy is inverted 0] to avoid this problem, then z/g is a few eV in mass and the 



constraint of neutrinoless double /3-decay that m^^ < 0.68 eV|TT|] cannot be satisfied. Given 
the manifold uncertainties of the r-process calculation in a hot-bubble scenario, we choose 
to disregard it in favor of the double /3-decay constraint. 

We view the large mixing but small mass splitting of z/^ and Vr as suggestive of their 
pseudo-Dirac origin. We also envisage a massless and a very light z/5 as two Majorana 
neutrinos with a small mixing. A unified seesaw model of these hierarchical masses would 



need a Dirac seesaw |12l] for the former pair with a minimum 4x4 matrix and a Majorana 



seesaw 



0] for the latter pair with a minimum 3x3 matrix. An additional small mixing 
between these two sectors is necessary for understanding the LSND results and can only 
come from both matrices being submatrices of one 7x7 neutrino mass matrix. 

It may appear difficult at first sight to construct a model of the lepton sector generating 
naturally a mass matrix with the above requirements from a symmetry. However, as shown 
below, a reasonably simple model does emerge if one supplements the standard SU{2)l x 
U{1)y electroweak gauge symmetry with a discrete symmetry ||14|| which might be the 
product of a more fundamental underlying theory. This discrete symmetry will be broken 
spontaneously resulting in the appearance of domain walls. However, it is now known |p!5| 



that higher-dimensional operators, induced at the Planck scale, can make such domain walls 
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collapse very quickly after formation so that we need not consider this as a potential problem. 

We take the Z5 elements to be a;~^, u!~^, 1, cu, and a;^ with uj^ — 1. Let the three lepton 
families of left-handed electroweak doublets be denoted by {ua, l, with a — e, /i, r. Let the 
three right-handed charged lepton singlets laR be accompanied by four right-handed neutrino 
singlets {vaRi^SR}- The Z5 transformations of the leptons bearing subscripts e,iJ,,T,S are 
chosen to be 1, o;""^, o;^, u;~^ respectively. The scalar sector is assumed to consist of two 
doublets $1 = {(pf, 05); ^2 = (02^, 0°) ^ complex singlet x° transforming as 1, a;~^, and 
cu respectively. The charged lepton mass matrix linking I^l and Ii^r is now of the form 

' a d' 

Mi^ e b , (1) 
0c 

where the diagonal entries a,b,c come from the nonzero vacuum expectation value of 0° and 
the off-diagonal entries d, e come from that of 02 • The zeros of this matrix are protected 
at tree level by the assumed discrete symmetry. As it stands, this mass matrix shows 
possible — tl but very little cl — /il niixing. There could be substantial cr — /ir mixing, 
but that is not observable as far as vector gauge interactions are concerned. 

Turning to the neutrino sector, we find the mass matrix spanning PeL, ^^Li ^tLi ^cR, ^hR, 
i^rR, and usR to be given by 
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m5 






(2) 



where mi,m2,ms come from (0?), m^^m-j from (02), 7714,7719 from (x°), m^,m^ from (x'^), 
and Ml , M2 are allowed mass terms even in the absence of symmetry breaking. The zeros 
are again protected at tree level by the assumed discrete symmetry. 



Note first that Ai^ has one zero mass eigenvalue, corresponding to an eigenstate which 
is mostly Ue as we will show. In the absence of symmetry breaking, all the m's are zero 
and we have only Mi and M2, corresponding to one massive Majorana fermion and one 
massive Dirac fermion respectively. These allowed masses can be very heavy and will act as 
the anchors for the seesaw mechanisms which generate the requisite small neutrino masses 
of our model. Note that the 4x4 submatrix spanning u^l, ^f^L, ^tLi and usr is indeed 
identically zero, even in the presence of symmetry breaking. 

Consider now the situation where only m2,m3 are made nonzero in addition to Mi,M2. 
Then we have a Dirac seesaw mass m2m^/M2 for and z/,-, whereas Ve and vs remain 
massless. Thus M.^ has a global vector U{1) symmetry as well as two chiral U{1) symmetries 
in this case. If the other m's are also made nonzero, then these symmetries are broken 
except of course for that corresponding to the zero mass eigenvalue noted before. Hence it 
is natural[|16] to assume that these other m's are much smaller in magnitude than 7712, 
which are in turn much smaller than Mi, M2. As it turns out, we also need to make the ratio 
rrii/m^ small because it corresponds to — i^s niixing for explaining solar neutrino data. 
To summarize, we assume 

|mi| << |m4,5,6,7,8,9| << \m2,3\ « |^1,2|- (3) 

For very large Mi, M2, the seesaw reduction of A^^ yields the following 4x4 mass matrix 
spanning usR, ^eL, i^f^L, and UrL- 

ml 17111714 m^rriT Tn^niQ m2Tn^ 

rriim^ ml mim^ 1 m^m^ fn^m^ 

M2 m2m^ m2m3 

m^m-j mim-j m^ rn^mQ m2m^ 

Taking the m2m^/M2 entries in the above mass matrix to be dominant, it is easily seen that 
a further seesaw reduction yields a 2 x 2 matrix spanning only vsr and UeL with elements 



M' 



1 

'Ml 
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given by the upper- left-corner submatrix proportional to 1/Mi. To find the eigenvalues A' 
of M-l which we denote by — A/M1M2, we write down its characteristic equation: 

= X{X^ - X'^M2{ml + mi + mj) 

— XlM^lmlml + mlml + m\m^ -\- m\m\) -\- 2MiM2mim6(m4m5 -|- 777.3777.7)] 
-|- MlM2[m\m\{m\ + ml) + (77737774 — 77757777)^(7772 + ^^e)]}- (^) 

Using the mass hierarchy assumed in Eq. (3), the four eigenvalues are easily obtained: 



a; = 0, 

~ Ml' 



777 2 7773 
Mo 



X' 



7777 1 / 777 2 7771 7773777g \ 

2M[ ^ 2M2 V ^3 ^ 7772 / ' 



77727773 



7777 



/ 77727775 ^ 7773777g\ 



M2 2Mi 2M2 V 7773 7772 y 

The corresponding mass eigenstates are then related to the interaction eigenstates by 



(6) 
(7) 

(8) 
(9) 
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I/V2 






Vr _ 









1/V^ 


1/v^ 







(10) 



In Eqs. (6) to (10) we have consistently retained only the leading terms. 

We now have our desired pattern of neutrino masses and mixing. We see that and u. 
are pseudo-Dirac partners with mass difference squared given by 

2r77yr7727773 



577734 ~ 



(11) 



M1M2 

The singlet neutrino Vs is mostly V2 which has a small mass, whereas z/g is mostly Vi which 
is massless. Their mixing is given by mi/m^. Furthermore, i/^ oscillates to i/g with mixing 
given by mQ/m2 as shown by Eq. (10). For illustration, let Mi = M2 = 100 TeV, 7772 = 10 
MeV, 7773 = 25 MeV, 7774 = 0.5 MeV, 7775 = 0.4 MeV, 7777 = 0.6 MeV, and 7771 = 20 keV; then 
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the common mass of and Ur is m2m^/M2 = 2.5 eV, the mass of z/5 is m\/Mi = 2.5 x 10^'^ 
eV, 5m|4 = 1.8 x 10~^ eV^, Ve — vs mixing is mi/m4 = 0.04, V(. — v^ mixing is m^jm2 = 0.04, 
yielding an effective sin^ 29 = 6.4 x 10~^ in either case as indicated by the solar-neutrino and 
LSND data. 



With four light neutrinos, the nucleosynthesis bound||T3 of < 3.3 is an important 



constraint. Although us is a singlet neutrino, it mixes with u^. and may contribute signifi- 
cantly to through oscillations. However, for the matter-enhanced nonadiabatic — vs 



oscillations which explain the solar data, this is not a problem ||18||. Note that if we had made 
mi >> m4 instead, then us would be mostly massless and the resulting scenario would be 
excluded by nucleosynthesis. By the same token, if we had tried to let z/^ oscillate into us 
to explain the atmospheric neutrino data, it would also be in conflict with nucleosynthesis. 
These requirements are sufficient to pin down uniquely our model of four light neutrinos 
provided we restrict ourselves to only two-oscillator effects]^ and the viewpoint that two 
almost degenerate neutrinos should be pseudo-Dirac. Of course, this forces us to have max- 



imal mixing in the — sector and we must disregard part of the Frejus data|jl9[, but if 
a smaller effective sin^ 1Q (say 0.7) is desired, then the underlying symmetry as well as the 
mechanism for generating Jv['^ become much more complicated]^]. As it is, all we need is a 
discrete Z5 symmetry and the seesaw mechanism. 

Although z/e — v-r mixing is very small as given by Eq. (10), it may instead come from 
the charged-lepton mass matrix of Eq. (1). Hence our model can also accommodate such an 
effect. Because two different Higgs doublets contribute to the lepton mass matrices, flavor- 
changing neutral-current processes do occur via the exchange of scalar bosons. For example, 
yU — > 67 is possible, but because all the Yukawa couplings are suppressed in this model, these 
are all negligible as far as present experimental bounds are concerned. 

Regarding the Higgs sector which consists of two doublets and a singlet, it can be shown 
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that the imposition of our discrete symmetry actually results in a continuous U(l) sym- 
metry which is of course broken in the Yukawa sector as x° couples both to lyeRi^SR and 
^tR^tR- This means that a pseudo-Goldstone boson appears with a mass of order m^/ATr 
which is too small to be compatible with present data. To avoid this problem, a simple 
solution is to enlarge the scalar sector with a real singlet rjo and a neutral complex singlet ri2 
transforming under as 1 and cu^ respectively. In order that they do not couple directly to 
the leptons, they are also assumed to be odd under an extra discrete Z2 symmetry. Because 
of the newly allowed terms rj2rj2Xi V0V2XXJ a-^id r/o^o, there is no longer any unwanted U(l) 
symmetry in this case and all the scalar bosons can be heavy. 

The 4x4 neutrino mass matrix given by Eq. (4) is similar but not identical to those 
of Ref . [7] . We assume that i/^ and Ut- to be pseudo-Dirac partners whereas Peltoniemi and 
Valle took i/^ and a linear combination of Uj- and Ug, while Caldwell and Mohapatra did 
not consider them as pseudo-Dirac partners at all. Our model also differs in that we have 
a massless eigenvalue in the neutrino mass matrix and they do not. However, the most 
important difference is in how we realize the desired form of the mass matrix. We use a 
simple discrete Z^ symmetry and the seesaw mechanism whereas both papers of Ref. [7] 
require more complicated symmetries as well as radiative mechanisms for mass generation. 
Hence their scalar sectors are much more involved and contain many particles of exotic 
hypercharge whereas we have only the usual doublets and neutral singlets. Presumably, these 
are experimentally accessible at the electroweak energy scale and would serve as discriminants 
of one model against another. 

In conclusion, we have shown that with all present desiderata, a simple model of four 
light neutrinos can be constructed with pairwise oscillations explaining the solar neutrino 
deficit {vg — vs), the atmospheric neutrino anomaly {v^ — Vr), and the LSND observation 
{I'l^ — i'e)- If the Sm^ of the last experiment is indeed 6 eV^, then there is also the cosmological 
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connection of neutrinos as candidates for hot dark matter, though that would be in marginal 
conflict with another neutrino experiment [pO|]. If the experimental inputs chosen by us do 
indeed stand the test of time, the matrix of Eq. (2) with the hierarchy of Eq. (3) could well 
be the harbinger of a complete theory of neutrino masses and mixing. 
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